APPROACH TO EQUILIBRIUM FOR THE STOCHASTIC NLS 



J. L. Lebowitz, Ph. Mounaix and W.-M. Wang 



Abstract. We study the approach to equilibrium, described by a Gibbs measure, for a system on a d- 
dimensional torus evolving according to a stochastic nonlinear Schrodinger equation (SNLS) with a high 
frequency truncation. We prove exponential approach to the truncated Gibbs measure both for the focusing 
and defocusing cases when the dynamics is constrained via suitable boundary conditions to regions of the 
Fourier space where the Hamiltonian is convex. Our method is based on establishing a spectral gap for the 
non self-adjoint Fokker-Planck operator governing the time evolution of the measure, which is uniform in 
the frequency truncation TV. The limit — > oo is discussed. 
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1. Introduction 



The nonlinear Sclirodinger equation is used to model a variety of phenomena in different physical 

systems, see e.g., [LRS, MCL] and references therein. It comes in two versions: the focusing and the 
defocusing cases. These correspond respectively to the sign of A in front of the nonlinear term in the 
Hamiltonian H of the (isolated) system, 

H:=H{u,u)= [ {\Vu\'^ +m\u\^)dx+— [ \u\Pdx, (1.1) 

where p is a positive even integer and m > 0. Here u is a complex valued field (function) on the d-torus: 
T'' := [0,277]"^ after identifying the end points. The time evolution of u satisfies the equation: 

.du dH 

^ dt ^ du (1.2) 
= —Au + mu + X\u\^~'^u. 

(Note that (1.1) differs from the H in [LRS, MCL] by a factor of 2.) 

In the defocusing case, A is positive, H is bounded from below and one then has global existence and 
uniqueness of smooth solutions in d < 4 for appropriate p; more precisely, for p < 4 in d = 4, p < 6 in 
d = 3 and all p in d < 2. This is the case of interest in constructive quantum field theory [GJ]. 

The situation is very different for the focusing case, A < 0. In this case H is not bounded from 

below and there exists a critical pc, Pc = 2 + 4/rf, such that for p < pc, there are unique global smooth 
solutions; for p = Pc there are unique global smooth solutions with small norm. For p > Pc but 
P < 2d/(d — 2), one needs to have small data to have unique global solutions, cf. e.g., [Bo3]. 

This dichotomy carries over to the possible existence of a (non-truncated) Gibbs measure with density 

M = exp [-^ / — : \u{x)\P : dx], (1.3) 
Jt<* P 

with respect to the free field Gaussian measure, where /3 > and : : indicates Wick ordering, which is 
needed for d > 2 and is a normalization factor. 

These measures exist for the defocusing case, A > for all p in d = 1 (without Wick ordering) and for 
p < 6 in d = 2, and p < 4 in d = 3 [GJ]. Furthermore, despite the fact that this measure is supported 
on rough functions, Bourgain showed that it is invariant under the dynamics given by (1.2) for d < 2 
[Bol, 2]. This means that the dynamics can be defined (after Wick ordering modification in d = 2) on 
a set of full measure with respect to /i. 

The focusing case is more delicate. Since H is unbounded from below, it is obvious that the measure 
/i cannot exist without some restrictions on its domain. It was shown by Lebowitz, Rose and Speer in 
[LRS] that in d = 1, the Gibbs measure exists for p = 4 when restricted to balls and that it exists for 
p = 6 with the additional condition of small |A|, cf. also [McKV, Z]. In d = 2, Jaffe [J] showed that the 
measure exists for p = 3 for real u when restricted to L'^ balls and after Wick ordering; while Brydges 
and Slade [BS] showed that this does not work when p = 4. 

In this paper, we shall deal with both the focusing and defocusing cases when the system is in contact 
with a heat reservoir of the Ornstein-Uhlenbeck type at reciprocal temperature /3. This problem was 
studied in [MCL], whose formulation we shall follow here, cf. also [LRS]. The time evolution is then 
described by the stochastic nonlinear Langevin equation: 



du+iiy + i){-A + m + X\u\P-'^)udt= y/W^dW, (1.4) 
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where > is the friction and is a complex Wiener process, i.e.. white noise in space and time. We 
note that for tlie moment, equation (1.4) is formal as it involves a very rough force on the right side. 

To give a meaning to (1.4), we work in Fourier space and restrict the number of Fourier modes n to 
be finite and constrain the Hamiltonian so that it remains bounded from below and is locally convex. 
We enforce this bound on the evolution by imposing appropriate Neumann boundary conditions in 
Fourier space. The problem then reduces to a finite dimensional Markov process on a compact space, 
whose stationary measure is given by a truncated Gibbs measure /ijv- This is formulated in sect. 2. We 
note that for fixed N, exponential approach to hn without any lower bound on the rate, follows from 
general probabilistic considerations, see [MCL]. 

The main new result of this paper is that the approach to is exponential with an exponent given 
by the size of the spectral gap for the non self-adjoint Fokker-Planck operator, which is uniform in the 
truncation N. This is proven in sects. 3-5. The proof is made possible by identifying the Fokker-Planck 
operator with a non self-adjoint Witten-Hodge Laplacian and a spectral analysis. 

We mention that a related Witten-Hodge Laplacian in the self-adjoint setting was introduced by 
Sjostrand [SI] to study correlation functions in equilibrium statistical mechanics and used by Wang [W] 
to study the parabolic Anderson model. 

The question of what happens to the truncated measure hn and the dynamics in the limit as ^ oo 
is discussed briefly in sect. 6. It appears that in d = 1, the limit will coincide with the given in (1.3), 
whenever the latter exists. The situation in d > 2 is less clear since Wick ordering takes away convexity. 

Acknowledgement. W.-M.W. thanks J. Sjostrand for several useful conversations. J.L. and W.-M.W. 
thank D. Brydges and G. Gallavotti for clarifications on the construction of the Gibbs measure. J.L. and 
Ph.M. thank P. Collet for many helpful discussions. J.L. is partially supported by the grants NSF DMR 
08021220 and AFOSR AF-FA9550-10; W.-M.W. is partially supported by the grant ANR-10-JCJC0109 
and the Qiushi professorial chair of mathematics at Zhejiang University. 

2. Truncated Fourier space formulation and the Theorem 

Let n e Z"^ be the dual variable to a; e T"^ and a = {a„} the dual of u: 

We rewrite equation (1.3) using Fourier series. Let $ be the corresponding sum in the Fourier space 
for H: 

<^ia)/(2TTf := $(a,a)/(27r)'^ 

= ^ (n2-hl)|a„|2-h^ X] aj,aj^...aj^_^aj^, ^^'-^^ 

where a„ e C, stands for |np and for simplicity we have set m = 1. We now restrict to u such that 
a„ = for |n| > N and make estimates uniform in N, cf. the remark above (5.1) in sect. 5. Since all 

estimates are uniform, we will generally omit TV as a subscript, except in the statements of the Theorem 
and its Corollary. For notational simplicity, we set (3 = 2 and rescale time by (27r)''. 

We have the following equation for a: 



(2.2) 



where the r„ are independent centered complex Gaussian white noises with 

(r„(t)f„(0) = 2<5(i-0- (2.3) 

The Fokker - Planck equation describes the time evolution of the probability density P{a, t) with 
initial distribution P{a, 0) for a process given by a stochastic PDE. Conjugating by e^* and setting 
P = e*P, the Fokker - Planck equation for (2.2) has the form: 

dtP{a,t) + LP{a,t) = 0, (2.4) 



where 



Here 



«2 

+i^||V$f -t/A$ + 7^j„ j/>0. (2.5) 



is the vector {da„^, da„^}, (2.6) 

i 

danda-n 



^^ = E^^' (2-7) 



and 



= 2Y,{da„^d-a^ - d-a^^daj (2.8) 



is the Hamiltonian vector field for the (finite dimensional) Hamiltonian $ and is anti self-adjoint. (All 

indices above are restricted to \n\ < N). Thus the Fokker - Planck operator L is non self-adjoint. This 
is a distinguishing feature of the present construction. 

In order for the Gibbs measure /in ~ e"^*"^"-* to be stationary for the Fokker-Planck equation (2.4) 
in the focusing case, we need to restrict to appropriate domains in a. Toward that end, we define the 
balls and the Hamiltonian H-ball of radius VB to be the set of a such that, respectively 

J2 («^' + l)l«n|' <B, B>0 (2.9) 

and 

|$(a)|/(27r)'^ <B, B > 0, (2.10) 

where $(a) is defined in (2.1). 

Since we are considering u such that a„ = for |n| > N, the sums in (2.9, 2.10) are restricted to 
\n\ < N. We take as our domain the Hamiltonian H-balls and impose adapted Neumann boundary 
conditions on dM - the energy surface: 

r/-(-!/V$-i/V + /i$)P|aH = 0, (2.11) 

where r] is the normal direction, assuming that it is well-defined and 

/l3. = {25a„$,-2aa„$} 

is a vector, displayed using the same convention as in (2.6). On occasions, we will also call the 
Hamiltonian vector field. It is important to remark that since is tangential, the above boundary 
condition is equivalent to 

^ . (_j.V$ - z/V)P|aH = 0, (2.12) 

which plays a pivotal role in establishing a spectral gap uniform in A''. 

In the focusing case, the H-ball is not necessarily a connected set. The theorem below addresses the 
component that contains the origin. We call it the Ho-ball. In the defocusing case, Hq coincides with 
H. 

Our main results are : 



Theorem. For d = 1, 2, p E 2N arbitrary, d = 3, p = 4, 6 and d = A, p = A, there exist C, c > such 
that uniformly in N, the Fokker-Planck semi-group satisfies 

||g-tLiv _ no||£(i,2(]H[(,),I,2(Ho)) < Ce'"*, t>0 

for sufficiently small \X\, in a given Hamiltonian Mo-balls with the Neumann boundary condition (2.11) 
on the Hamiltonian surface OMq, where Ho is the projection onto the unique ground state of L^: 

and Cn is such that the truncated Gibbs measure fiN '■= V'ir normalized: 

I^N Yl dandcLn = 1, 
\n\<N 

and we have indicated explicitly the truncation N on the i>(o) defined in (2.1). 

Corollary. In the defocusing case, the Theorem holds for all X > under the same conditions on d 
and p. 

We note that in the Id focusing case, for p = 4 and 6, the Theorem remains vahd with the Hamiltonian 
Ho-ball replaced by an L^-ball, cf. sect. 6. This is because the Hamiltonian vector field is also 
tangential to the L^-sphere and (2.12) remains valid on the sphere. Without dissipation and forcing, 
this Hamiltonian geometry leads to the conservation of mass, i.e., the norm and energy, essential 
for the global existence and uniqueness of solutions for subcritical (or critical) dispersive NLS. With 
dissipation and forcing, the conservation laws arc no longer there, but it is important to recognize that 
the geometry remains and is precisely what enables the present construction. The Id focusing case has 
particular interest as the limiting measure fj. is known to exist [LRS], cf. sect. 6. 

Previous related results on stochastic NLS, e.g. [EH, KuS, O] pertain either to the focusing L^- 
subcritical case, namely p < 4, the Id focusing case, where there are no blow-up solutions, or to the 
defocusing case. The forces there are smooth in the spatial variable x. So one works directly with 
the (non-truncated) SNLS. Generally speaking, this approach does not yield the explicit form of the 
invariant measure. For results in the purely dissipative case, cf. [DD]. 

Equation (1.3) balances dissipation with a rough force. The main novelty of the present construction 

is to work directly in the phase - space (the cotangent bundle) using the Fokkcr - Planck operator 
and impose boundary conditions there. This enables us to smooth the rough force by convexity of the 
Hamiltonian and treat the critical cases. 

3. Sobolev embedding and convex Hamiltonian 

We establish convexity of $ uniformly in N and prove 

Proposition. Under the same conditions on d and p as in the Theorem, the Hessian of the truncated 
$ satisfies 

cl < Hess $ := < CI, (3.1) 

in the ball defined in (2.9), where c, C > are constants which do not depend on N, if |A| is 
sufficiently small. Here Hess is considered as an operator on i'^{[—N,N]'^) x £'^{[—N,N]'^). If ^ is 
defocusing, then A can be taken to be 1 and C = C{B). 
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Proof. We write the matrix operator in the form = Mi + M2, where 



[[da AM J 



(3-2) 

\[[daAM 

and [[ ]] denotes the matrix of second order partial derivatives. 
Now to prove (3.1), it suffices to prove that 

{v,^"v)>c{v,v) (3.3) 

for all V = {wj,Wj} e x t^. The quantity {v,^"v) is a sum of terms of the form Miji~WjWk or 
M2.jkWjWk and their complex conjugates. Let W be the function with Fourier coefficients Wj and u 
the function with Fourier coefficients aj (|j| < A'^). We observe that since $ is a sum of homogeneous 
polynomials in aj,ak, the terms that occur in the sum are exactly the same as the terms in the evaluation 
of the following integrals: 

[m%^ := [\vw[\h + [mh = + (3-4) 

\\u^p/'^-'^^W[\l2, J W'^u^[u[^P-'^Ux and its complex conjugate, (3.5) 

where here and below J ■ dx denotes integration over tlie torus T"^. 
More precisely, 

{v, ^"v) = 2\\W\\]j^ + \p\\u^P'^-^^W\\l2 + X{p - 2)Re J W\^[u[P-^dx. (3.6) 

Using now the inequalities 

< ||^(p/2-i)^||2 ^ /J _ 2/p)Re [ W^u'^[u[P-^dx < 2\\u^P/^-^^W\\l2, (3.7) 
P J 

we have in the focusing case, 

2\\W\\%^ - 2\X[p\\u^P/^-^^W\\h < {v, ^"v) < 2\\W\\%^ + 2[X[p\\u^P/^-^^W\\l2 (3.8) 
and in the defocusing case 

2[[W[\l. < (v, $"^;) < 2[\W[\l. + 2[X[p[\u(P/'-'^W[\l2. (3.9) 
Standard Sobolev embedding on T'', namely 

II • IIlp < C'd.pll • ||iir<i/2(i-2/p) (3.10) 
then gives under the same conditions as in the Theorem on d and p that 

||«(^/2-i)ty||i. < C,J\u\\P„-,'\\W\\h. (3.11) 

Using (3.11) in (3.8, 3.9) and also smallness of |A| in (3.8), proves the proposition. □ 



4. The Witten-Hodge Laplacicin and spectral gap 

We now prove that the Fokker-Planck operator L in (2.5) restricted to Hamiltonian H-balls defined 
in (2.9) has a spectral gap uniform in N when the Hamihonian is convex. Wc use the Wittcn-Hodge 
Laplacian formulation introduced in a related context in [SI]. We summarize some of the basic notions 
in the Appendix using the self-adjoint setting. Here we show that it corresponds to the operator L and 
gives the desired spectral gap. 

The non self-adjoint Fokker-Planck operator 
Let d be the exterior differentiation 

d:=Y,dbJb^, 

n 

where 6„ stands for a„ and a„ and 

:= e-*rfe* = ^idb„ + 4„$)^^&^ (4.1) 

n 

We note that here we work on with A = [-N, iV]<^, which could be identified with (M2)A M?W. 
Let 

where the notation is the same as in (3.2) with each bloc being one of the four possible sectors: 
uu, uu, uu, uu. Define the (formal) adjoint of with respect to A to be: 

n 

(Properly speaking, A is a map from the cotangent space (C^"-^'-^] ) to the tangent space 
The Witten-Hodge Laplacian is then defined as 

— A$_A = d*^^dis, -h 

cf. the paper of Bismut [Bi] for the general construction and also the Appendix in sect. 7. 

The reason that wc introduce the Wittcn-Hodge Laplacian is that when restricting to scalar functions, 
the 0-forms, it is precisely the Fokker - Planck operator L in (2.5) and we have the following identities: 

-A^A = di^d^ =dloAod^ = L, (4.2) 

where 

n 

We will also need it on 1-forms: 

-A^^_)^ = -Ag)^ (g) I + 2$" o A\ (4.3) 

cf., [HHS, S2]. 

Remark. This is occasionally dubbed the supersymmetric approach, cf. e.g. [TT-NK] for the physics 
literature. 

Let the Hamiltonian H-ball be as defined in (2.10) and as before, Ho-ball the connected component 
containing the origin. The following basic spectral characterization of the Fokker -Planck operator 
follows readily from (4.2). 
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Lemma 1. Assume Mo is a convex set. The spectrum of L in this ball with the Neumann boundary 
condition (2.11) is contained in the sector 

{k e C\ReK > u\Im,K\}, (4.4) 

and is a simple eigenvalue. 

Proof. The spectrum of L, a{L) is contained in the numerical range of L: 

{{u, Lu) e C|u in the domain of L}. 
Prom (4.2) for all u in the domain of L satisfying the Neumann boundary condition (2.11), 

{u, Lu) = {u, — lS.'^^\u) = V i \dis,u\'^ dada + i Im. i (d^u)^dada, (4-5) 

where {d<s>u)'^ denotes the sum of the square of components, with each component defined as in (4.1). 
Equation (4.5) implies (4.4), since (4.5) equals zero if and only if dg,u = 0, this means u is a multiple 
of e~*, which concludes the proof. □ 

Using the Witten-Hodge Laplacian (4.3) on 1-forms, we have, moreover, the following spectral gap 
lemma, essential for the proof of the Theorem. 

Lemma 2. Assume that in the Mo-ball, $" satisfies cl < $" < CI for some c,C > 0. Then for som,e 
Co = co(c) > 0, the Fokker - Planck operator L with Neumann boundary conditions has the following 
properties uniformly in N: 

(i) {o-(L)\{0}} n e C\Re z<co} = 0; 

(ii) (up to constants) e~* is the unique ground state with eigenvalue 0. 
Let Ho be the projection onto the normalized ground state. Then 

(Hi) \\{L — 2;)u|| > (co — Re z)\\u\^ if u = {1 — no)u. 

Proof. Assume u is an eigenfunction of the Fokker -Planck operator L with Neumann boundary condi- 
tions with eigenvalue k^O, i.e., 

Lu = -A^°^^M = Ku, 0. (4.6) 

Then from Lemma 1, d$w ^ 0. Operating on equation (4.6) using d$ and taking the scalar product 
with d^u , we have 

{d^u, [d<s>dy^]d^u) = K{d<!>u, d<!>u) . 
Writing v for d^u, we have equivalently 

(v, -A^^t;) = {v, [-Ai°)^ ® I + 2$" o A» = Kiv, v), (4.7) 

where we used (4.2). 

Taking the real part of (4.6), we obtain 

Re K\\vf = {v, [-A^°^ O l]v) + 2u{v, ^"v), 

where — A^^ = d'^d^ is the self-adjoint Laplacian on 0-forms. 



Using the convexity of we then obtain (i) with cq = 2cv > uniformly in N. Here we also used 
the fact that the self-adjoint Laplacian — A^^ has the same Neumann boundary condition as — A^^^. 
We remark that in fact stronger results are known under appropriate conditions, namely 

a(-Ai%)\{0} C a(-A«^), 

cf. [HHS]. 

(ii) follows from (i) and Lemma 1. To prove (iii), we write 

||(L-^H||H| = ||(-Ai%-^H||H| > \{-A'i^^-z)u,u)\ 

>\Re{{-A<iJ^-z)u,u)\ ^^g^ 

= |((-Ai°)- Re^K«)| 
> (co — Re z){u,u) 



for all u = (1 — Ilo)u. 

Here we used self-adjointness of — A^\ the fact that up to constants e~* is also the unique ground 
state of — A^\ i.e., 

no(-Af) = no(-AWj 

and eigenfunction (of — A^') expansion of u. □ 

Remark. Since L is non self-adjoint, it is no longer true that the resolvent at z is bounded above by the 
inverse of the distance of z to the spectrum. This is because of the non-commutativity of the self-adjoint 
and anti self-adjoint components and hence the necessity of the type of arguments in (4.8). 



5. Exponential approach to equilibrium 

Proof of the Theorem. We only need to verify the conditions in Lemmas 1 and 2, namely convexity of 
the Hamiltonian in Hp. The rest will follow by applying these two lemmas and contour integration. 
Recall that the H'^ and the Hamiltonian H-balls arc as defined in (2.9, 2.10) and we restrict to functions 
u such that its Fourier coefficients a„ = for |n| > N. So the sums below are to be understood as: 

E-E= E; 

nGZ"* \n\<N 

in other words, the B below is fixed independent of N. Therefore the radius of the various "balls" are 
fixed. It is only the dimension of the balls that varies with N. 

Remark. For most of the applications, u = ttu, where tt is the projection onto the first 'W Fourier 
modes while u has full Fourier support and is in at least one of the "balls" . 

Define the Hamiltonian H-ball as before: 

2A 

|$(a)/(27r)'^| =1 Y^{n^ + l)|a„|2 + — ^ aj,aj,...aj^_,ajj 

^ ji-j2...+jp-i-jp=o (5.1) 

< B. 

Using the results of sect. 3, we have that if a is in the above H-ball, then for the defocusing case, A > 0, 
$ is convex and {v, ^"v) < C. So H is a convex set with the well-defined boundary 



dU = {ama)/{2nf = B}, 
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and 

^{n^ + l)\an\^ <B 

for a in H. 

For the focusing case, A < 0, assume we look at the connected component containing the origin, 
namely Hq 9 0. Let A be the set 

^ = {a|^(n2 + l)|a„|2<5B}, (5.2) 

and 

Ma := Ho n A. 

The Proposition then gives that $ is convex in A and moreover $(•) is equivalent to the norm: 



< (1 - c\x\BnmH, < V §^ < + c\x\Bn\\f\\H, < 2|i/ik„ (5.3) 

for |A| <C 1 and where / is the function with Fourier coefficients a G A. 
Equation (5.3) together with (5.2) give that Ho is strictly contained in A: 

Ma = Ho. 

So Ho is a convex set with the well-defined boundary 

dMo = {a\ $(o)/(27r)'^ = B} c A. 

We now evaluate the semi-group e~*^ using the contour F = Fi +F2, where Fi is compact enclosing 
and Re Fi < c' < Co, the co in Lemma 2; and F2 is defined by 

{z\ |Im z\ = |Re z\'^ + D, if Re ^ > c' and |Im z\ < c'^ + D if Re z = c' for some D > 0}, 



so that 



= — { e-*'(z- ly^dz 



= — i e-*'{z- L)-'^dz+ — { e-^'^iz- L)-^dz. 

It follows from Lemmas 1 and 2 that the above contours lie in the resolvent set and the integrals are 
well-defined, (i, ii) of Lemma 2 proceed to give that the first term is Ho and (iii) gives the exponential 
estimate for the second term, uniformly in N and concludes the proof. □ 

Proof of Corollary. This follows from global convexity of which is moreover uniform in the radius 
^/B of the H-ball, since ball is equivalent to H-ball in the defocusing case. Since 

^ Obn 



where 6„ = a„ and a„, let u be the function with Fourier coefficients {a„}, we have that 



where we used convexity for the lower bound, (3.10) and the restrictions on d and p. So the energy 
surfaces 9H is well-defined if u € fl H^. Using the density of in H^, we reach the conclusion of 
the Corollary. □ 
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6. The limit N ^ oo 



As discussed in tlie introduction, the existence of a Gibbs measure n ~ exp[— for the H in (2.1) 
is a problem which has been stiidicd extensively. It corresponds in = 1 to the problem of the behavior 
of the limit A'' — > oo of the measure /j-n ~ exp[— /3$jv]- These measures exist for all p in the defocusing 
case and for p < 6 for the focusing case under some restrictions. In d > 2, even the defocusing case 
requires Wick ordering, see (1.3), and it does not exist for the focusing case. 

In d = 2, the Hamiltonian dynamics corresponding to = 0, have been modified by Bourgain to 
include Wick ordering and shown, as mentioned earlier, to be well defined for the defocusing case when 
p = 4. Whether one can make a similar modification to the Langevin dynamics, 2/ > 0, is an open and 
intriguing question. We note however that the Wick ordered nonlinear term : : is not convex. So 
the argument in this paper leading to the existence of a spectral gap would not hold for the modified 
dynamics. Of course the result might still be true. But this is another question. 

The question of what happens to the dynamics generated by (2.2) or the Fokker - Planck equation 
(2.4) in the limit A'^ — >^ oo is therefore of relevance here primarily in the case d = 1. We consider the 
defocusing and the focusing cases below. 

The defocusing SNLS 

In this case the estimates in the Corollary are not only uniform in N, but also in the radius \/B of 

the Hamiltonian ball, which is equivalent to the norm. For d = 1, letting B ^ oo and subsequently 
N ^ oo, should therefore lead to the limiting measure with density: 



with respect to the free field Gaussian measure. A full proof is under investigation. 
The Id focusing SNLS 

In the Id focusing case, the limiting measures when p = A and for small |A| also p — 6 are known 
to exist [LRS] and are invariant [Bol] under the corresponding dispersive NLS dynamics, i.e., (1.3) 
when = 0. When v > 0, the Theorem indicates a spectral gap uniform in A^, albeit restricted to a 
Hamiltonian ball in which it is convex. 

For fixed A^, this restriction to a Hamiltonian ball can be replaced by a restriction to an ball with 
Neumann boundary conditions, as mentioned after the Theorem in sect. 2, since the Hamiltonian vector 
field is also tangential to the sphere. We then have that for both p = 4 and 6, the Hamiltonian is 
convcix (uniformly in A^) for small |A| depending on the radius of the ball using Sobolev embedding 
and interpolation. In the focusing case, the Hessian of the Hamiltonian is bounded above by a constant 
(which only depends on N). So using the same arguments as in sects. 4 and 5, the Fokker-Planck 
operator has a spectral gap uniform in N. 

This Neumann restriction is a priori compatible with the ensuing N ^ oo limit as the Brownian 
motion in d = 1 is supported in H^^^ C and should lead to the limiting measure constructed in 
[LRS]. The question would then concern its invariance with respect to the SNLS dynamics in (1.3). 
The general procedure could be akin to that in [Bol] on the corresponding Hamiltonian dynamics. 
Alternatively one could try to consider the limit A" — >■ oo of the equation (2.4). See [H] for some results 
in a related but finite N setting. 



Below we make a short introduction to the Witten-Hodge complex (originally introduced in [Wi]) 
using the self-adjoint setting. This is because the structure remains the same and the notations are a 
bit simpler, cf. Chap. 11, in particular Chap. 11.4 in [CFKS]. 
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7. Appendix: The Witten-Hodge complex 



Toward that purpose, let </> e C°°(IR^;IR). Let d be the usual exterior differentiation: 

N 



and 

d^ = e-'I'de'^ = d + d(t)^ = J2 '^^i ' C^-l) 

where 

(For the calculus of differential forms, see for example [Sp].) 

If / is a form of degree m, then d^/ is a form of degree m + 1. For example, if / is a 0-form, i.e., a 
scalar function in C°°(R^;R), then 

N 

is a 1-form, which we may identify with a vector valued function F with the components: 

Fj{x) = {zjf){x), 

i.e., a function in C°°(R^; K^). (We note that when = 0, d^f = df, which is just the usual differential 
of /.) If / is a 1-form, / = J2j fjdxj then 



d,j,f = '^{zifj)dxi A dxj 



%<3 

is a 2-form, which we may identify with an N x N antisymmetric matrix function M with the entries: 

M,j{x) = -Mji{x) = (z,/,)(.t), 

I.e., a function in C°°(K^; A M^). The operator Zj can be seen as an annihilation operator in view 
of its action on e~*, namely 

zje-'^ = 0, for j = 1,2,..., AT. 

On the space of m-forms, one may define an scalar product for two m-forms uj and v: (w,u). 
Specializing to v = d^f, where / is an (m — 1)- form, we define the adjoint operator of d^ so that: 

{w,d^f) = {dlujj). 

This gives the formal adjoint to be 

TV 

d; = e'l>d*e-'^ = J2^jdx\, (7.3) 

where 

■' dxj dxj 

and J is the usual contraction which lowers the degree of the forms. 
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If / is a form of degree m, then d^/ is a form of degree m — 1. For example, if 

N 

is a 1-form, then 

AT 

is a 0-form, i.e., a scalar function in C°°(]R^;]R). (We note that when ^ = 0, d^f = d*f is just the 
usual divergence of /.) If / is a 0-form, then d'^f = 0. 

The operator z* can be viewed as a creation operator. For example, when N = 1 and (f) = x^, 
z* := Zj generates the first Hermite polynomial. We have the commutation relation: 

[zj,zl]=2djdkcP, (7.5) 

which plays an important role. We check easily that indeed 

d^d^ = d%dl = 0. (7.6) 

Using dcj,, d*^, we define the Witten Laplacian, 

-A^ = + rf^rf;, (7.7) 



on C°°(R^; A%^) {1<1<N), where denotes the anti-symmetric tensor product of with 

AM^ := M^, A^M^ := AR^ etc. For example, if £ = 2, then an element of C~(M^; AM^) is an 
N X N anti-symmetric matrix valued function mentioned earlier and when £ = N,it is the determinant. 

Notice that 

d^A^ = A^d^, d^A^ = A^dl (7.8) 
by using (7.6). If we let A^^' be the restriction of A^ to forms of degree £, we obtain more precisely: 

d,A^:' = A^r^d,, d;Af ^) = Afd;. (7.9) 

Wc remark that the above construction is similar to that of Hodge Laplacian which corresponds 
to taking ^ = 0. (For a quick overview of the analytical aspects of Hodge theory, see Chap 11.3 in 
[CFKS]. ) We have explicitly 

-Af = dld^, = J2 ^Pj = - E + ll^-^ll' - Hess <l>. (7.10) 

j=l j=l 3 

For example, if is a non-dcgcncratc quadratic form, then — A^"-* is an A'^-dimcnsional harmonic os- 
cillator. Zj and z* are just the annihilation and creation operators for the harmonic oscillator. More 
generally, we have 

-A^ = E E ^j^dXjdxl + E E ^^jdxldx^ 

= E E ^ZjidXjdxl + dx[dx'>) + [zj, z^]dx^dx{ 

= E + ^ E E^^^^^^^i^^^^^i ^4 

= -Af 01 + 2^ J2(9,A>^,j,)dx;>dxl 
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where to obtain the third hne from the second, we used (7.5). In particular, with the identification of 
1-forms with valued functions, we obtain 



-AJ,'' =-Af ®I + 2.^". (7.12) 

Since formally (— A^^-*?*, u) > 0, under appropriate conditions on (f> at infinity, we can define — A^^-* 
as a self-adjoint operator and that — A^^'' has compact resolvent, cf. [SI]. Moreover — A^^-* has discrete 
spectrum contained in [0, oo). The lowest eigenvalue of — A^^-* is zero and a corresponding eigenfuction is 
e~'^, since this function is annihilated by rf^. This eigenvalue is simple, for if u is another eigenfunction 
associated to the same eigenvalue, then = {—A^^^u,u) = \\d^u\\^ and hence d^u = 0, which means 
precisely that m is a multiple of e~'^. 

Using (7.8), we obtain the following intertwining property of the spectra: 

a(-Aj^)\{0}cc7(-A«). (7.13) 
This is because if / is an eigenfunction of — A^°^ 

-Aff = Kf (7.14) 
with K > 0, then operating on (7.14) with d^, we obtain 

-d^Aff = -{d^d%)d^f = -Af{d^f) = Kd^f. 

So if k; ^ 0, then d^f ^ is an eigenform for —A^^\ which is the statement in (7.13). 

This is in fact the main reason that we introduced — A^"^' . Using (7. 12) , we then obtain that cr(— A^'^' ) 
has a spectral gap if </> is strictly convex. We end this self-adjoint introduction by remarking that if e~'^ 
is the eigenfunction of the Schrodinger operator 
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for the lowest eigenvalue k, then 



Prom this point of view, the Witten-Hodge Laplacian can be seen as a natural generalization of harmonic 
oscillators when cf) is quadratic. 
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